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Abstract. We study positional coincidences between 
gamma-ray bursts in the BATSE 3B catalogue and 
quasars/AGN taken from the Veron-Cetty & Veron com- 
pilation. For most classes of AGN, for BL Lac objects 
and for radio-loud quasars, we find no excess coincidences 
above random expectation, and we give upper limits for 
their burst emission rate. 

However, surprising evidence is found for a positional 
correlation between gamma-ray bursts and radio-quiet 
quasars. A total of 134 selected bursts with a position 
error radius <1?8 and the nearest and intrinsically bright- 
est radio-quiet quasars (M v < —24.2 and z < 1.0) have a 
probability of >99.7% to be associated with each other. 
An analysis of a smaller sample of well- localized interplan- 
etary network gamma-ray burst positions supports this 
result. This correlation strongly favours the cosmological 
origin of gamma-ray bursts and enables to estimate its 
distance scale. The average luminosity of those gamma- 
ray bursts which we associate directly with radio-quiet 
quasars is of the order of 4 x 10 52 erg (assuming isotropic 
emission) . 

Key words: Gamma-rays: bursts, AGN: general, 
Quasars: general 



1. Introduction 

The nature of gamma-ray bursts (GRBs) remains an 
enigma, even nearly 30 years after their discovery with 
the Vela satellites in 1967 (Klebesadel et al. 1973). The 
bursts detected with the Burst and Transient Source Ex- 
periment (BATSE) on board the Compton Gamma-Ray 
Observatory (CGRO) show that their positions are com- 
pletely consistent with isotropy. For the subsample of 1005 
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GRBs in the 3B BATSE catalogue the Galactic dipole mo- 
ment is within 0.9 a of the value expected for isotropy, and 
the observed Galactic quadrupole moment deviates only 
0.3 a (Briggs et al. 1996) from isotropy. 

Thus, GRBs are distributed significantly more 
isotropic than any known Galactic population and most of 
the Galactic models are no longer in agreement with the 
data (cf. e.g. Briggs et al. 1996, Hartmann et al. 1995). 
A spherical dark matter halo model could still be con- 
sistent with the measured burst positions, but requires a 
core radius larger than necessary to explain the Galaxy's 
rotation curve (Hakkila et al. 1994; Caldwell & Ostriker 
1981; Briggs et al. 1996). Extended halo models have to 
place the typical BATSE GRB farther away than about 
100 kpc in order to obtain the measured isotropy. On the 
other hand, the observed bursts must be closer than about 
350 kpc, because no excess of bursts is observed in the di- 
rection of M31. Neutron stars born with velocities exceed- 
ing the local Galactic escape velocity have been proposed 
to match these requirements (Podsiadlowski etal. 1995), 
but in order to obtain isotropy of bursts, only a very small 
fraction of low- velocity neutron stars are allowed to burst 
in this model. Moreover, a delayed turn-on is necessary to 
avoid a concentration of bursts to the Galactic disk (see 
e.g. Li & Dernier 1992; Briggs et al. 1996). 

The alternative possibility, a burst location at cos- 
mological distances, does not require to postulate addi- 
tional restrictions (Paczyhski 1986). Also, the observed 
inhomogeneity of the bursts' radial distribution can be 
regarded as supportive evidence for this idea. The num- 
ber of faint bursts shows a deficit with respect to the 
power law of slope —3/2 of a homogeneous distribution 
in Euclidean space (Schmidt et al. 1988). This can be ex- 
plained naturally if faint bursts are located at high redshift 
(Paczyhski 1991; Dermer 1992; Piran 1992). There have 
been many attempts to estimate the redshift range of GRB 
sources from their brightness distribution (e.g. Fenimore 
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etal. 1993b; Fenimore & Blom 1995). Whereas a redshift 
of about 0.1 for the nearest bursts seems likely, only a 
very crude estimate of the redshift of the faintest bursts 
can be made (z = 0.8 through z = 3.0). Apart from the 
uncertainty in the cosmological parameters, such studies 
suffer from the unknown evolution function of the burst 
sources and the corresponding emission rate (see Sect. 4.1) 
as function of redshift (Rutledge et al. 1995). 

The realization that gamma-ray measurements alone 
may be not sufficient to solve the GRB mystery has 
prompted increasing attempts for counterpart searches of 
GRBs at various wavelengths. Two different methods are 
possible to search for GRB counterparts: 

(1) Because of the typically large errors of burst posi- 
tions and the huge number of quiescent objects within 
each GRB error box, one method consists of simultaneous 
observations at other wavelengths (X-ray, optical, IR, ra- 
dio) to search for transient sources during the GRB event 
(see recent reviews by Hartmann 1995, Greiner 1995). As a 
result of these studies, there is solid evidence that gamma- 
ray bursts do not produce simultaneous optical flashes 
brighter than about 5 th mag for an assumed duration of 
1 sec (Greiner et al. 1996) or simultaneous radio emission 
stronger than 10~ 13 erg cm~ 2 s _1 (Cortiglioni ctal. 1981; 
Inzani etal. 1982). 

(2) The second, more statistical approach is the search 
for positional coincidences of GRB locations with com- 
parison samples of any other known class of astronomical 
object. Various types of catalogued objects, both Galac- 
tic and extragalactic, were compared with the BATSE 
gamma-ray error boxes, but any correlation between these 
and the bursts was compatible with random expectation, 
(i) Webber et al. (1995) looked for catalogued counter- 
parts of 60 bursts with the smallest error boxes available 
(~ 0.25 dcg 2 ) and found no position coincidences above 
random expectation for extragalactic objects (see, how- 
ever, Sect. 5.2.2. for our investigation of the latter sam- 
ple), (ii) Gorosabel et al. (1995) found no significant corre- 
lation between 44 GRBs localized by the WATCH instru- 
ment to less than one square degree and objects (galac- 
tic and extragalactic) listed in 33 different catalogues, 
(iii) Extremely rare events in external galaxies, for ex- 
ample merging neutron stars, are one possibility to ex- 
plain an extragalactic origin of GRBs (Eichlcr etal. 1989; 
Paczyhski 1991). Under such a hypothesis a correlation 
between burst positions and the large-scale structure of 
matter should be expected. The concentration of nearby 
galaxies towards the super-galactic plane offers the possi- 
bility to search for counterparts up to a redshift of z ~ 0.1. 
Hartmann et al. (1996) found no compelling evidence for a 
concentration of GRBs to the super-galactic plane in the 
3B catalogue, (iv) With a likelihood method that com- 
pares the counts-in-cell distribution of gamma-ray bursts 
with the one expected from the known large-scale struc- 
ture of the Universe, Quashnock (1996) argued that the 
nearest bursts originate at a redshift z ~ 0.03 and the 



farthest ones at z > 0.25. (v) However, searches for a cor- 
relation between GRB and the direction of rich clusters 
of galaxies (situated right in the above redshift range) did 
not arrive at consistent results (Howard et al. 1993, Ne- 
miroff ctal. 1994, Deng et al. 1995, Kolatt & Piran 1996, 
Marani et al. 1996). From a comparison of the angular cor- 
relation functions of GRBs and that of extragalactic radio 
sources, Kurt & Romberg (1995) concluded that gamma- 
ray bursts do not belong to a normal population of galaxies 
with z < 1.5. 

If we assume that GRBs are emitted at cosmological 
distances and keep in mind that no clear correlation was 
found with either clusters of galaxies or the super-galactic 
plane then it seems unlikely that any normal constituents 
of galaxies can cause gamma-ray bursts. The most promis- 
ing sources with a different nature are quasars and active 
galactic nuclei (AGN). They had been suggested as GRB 
sources by Prilutski & Usov (1975) shortly after the dis- 
covery of GRBs (Klcbcsadel etal. 1973). 

Both, the number of bursts with relatively small er- 
ror boxes, and the number of identified AGN and quasars 
are increasing permanently Webber et al. (1995) used 
the 1989 version of the Veron-Cetty & Veron compilation 
of QSOs and AGN. The fact that the number of known 
quasars and AGN has doubled since then has motivated us 
to check their results with improved statistics. A further 
chance to find correlations arises if only a certain subclass 
of the AGN or quasars are counterparts of GRB. In that 
case the mean distance between the members of the sam- 
ple will increase and the error of the gamma-ray burst 
position will play a less important role. We have therefore 
not only repeated the work of Webber etal. (1995) with 
larger samples, but extended the search for correlations to 
various subclasses of AGN. 

Our paper is arranged as follows. In the first part of 
Sect. 2 the GRB sample and the methods used to search 
for correlations are introduced. In the second part we pro- 
vide the samples of AGNs and quasars as well as different 
subclasses. Our search for correlations between GRBs and 
various subclasses is given in the first part of Sect. 3. In 
the second part we study in more detail those subclasses 
which show a positive correlation. The final results are 
described in Sect. 4 and we elaborate on the one posi- 
tive correlation found. Finally, we discuss our results in 
the light of previous studies in Sect. 5 and conclude with 
some remarks in Sect. 6. 

2. The Samples and the Analysis Method 

2.1. The sample of observed bursts 

The 3B catalogue of GRBs (Meegan et al. 1996) provides 
the position of 1122 events, together with a statistical er- 
ror, r st , which is defined by the formal covariance matrix 
from a x 2 -test on the assumption of normal errors. The 
circle with this radius centered on the fitted position of 
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the burst corresponds to the 68% confidence ellipse. This 
error is based on the Poisson uncertainty in the BATSE 
measurement of the burst flux by each Large Area Detec- 
tor and is believed to have a Gaussian distribution (cf. 
the description of the 3B catalogue). For many bursts the 
additional systematic error of 1?6 is larger than the for- 
mal error. Its distribution is not very well known, but may 
have a more extended tail than a Gaussian (Pendleton et 
al. 1996). Combination of both errors results in a total 
error, R t , of 

Rt = v^) 2 + (1?6) 2 (1) 

which is an estimate of the 68% confidence interval for the 
burst location. 

In order to test a given source sample of QSO or 
AGN {"comparison sample" in what follows) for correla- 
tion with GRBs we calculate the number of coincidences. 
A comparison with the number of correlations obtained 
for simulated burst positions allows us to estimate the 
significance of the correlations obtained for the real burst 
positions. To test for any bias introduced by either the 
distribution of GRB positions or the non-homogeneous 
distribution of the comparison sample we will use three 
different methods. They differ in the way of calculating 
the number of coincidences and in the way the simulated 
burst positions are generated and are described in the next 
paragraphs. 

Due to the uneven exposure of the sky the burst distri- 
bution is not isotropic. However, since the exposure map 
is primarily symmetric with respect to the ecliptic (except 
the south atlantic anomaly) , and the quasar distribution is 
asymmetric with respect to the galactic plane (due to the 
search techniques), the affect on the coincidence statistics 
is negligible. 

2.2. Randomly simulated GRB samples 

We calculate the number of GRBs containing at least one 
object of a given comparison sample within their 68% con- 
fidence interval. To reduce the noise all bursts which have 
a positional distance d < Rt to a source in the comparison 
sample are counted only once. 

Until now there are no reports on the presence of any 
large-scale structure (like e.g. dipole or higher moments) 
in the GRB distribution. Therefore it appears reasonable 
to use randomly simulated control samples of GRBs to 
obtain the number of chance coincidences, and their dis- 
tribution. However, we should keep in mind that there 
is no ultimate proof for the assumption that randomly 
distributed positions reproduce the characteristics of the 
observed GRB distribution. To estimate the significance 
of the correlation we simulate 1000 different sets of ran- 
domly distributed GRBs and determine the corresponding 
number of coincidences, C. The number of coincidences 
expected by chance, < Cb >, is given by the mean over 
the 1000 values of C. 
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Fig. 1. The expected number of position coincidences between 
Gamma-ray bursts and sources are plotted against the num- 
ber of sources in the studied sample. The used burst sample is 
given through a maximal position uncertainty of Rt < 1?8 and 
contains 136 gamma-ray bursts. If GRBs and sources of the 
comparison sample were distributed isotropically on the sky, 
the expected number of position coincidences would asymp- 
totically approach 136 at more than «4000 sources. However, 
quasars in the Veron-Cetty & Veron catalog are distributed 
very irregularly with nearly no occurrences in the galactic plane 
while clustering in fields where deep searches have been per- 
formed (especially for high-redshift quasars). 

Further, the one thousand values of C represent the 
distribution for coincidences between the comparison sam- 
ple and simulated burst sample which (by definition) do 
not show any correlation (null- hypothesis). This distribu- 
tion will be used to apply statistics on the coincidences 
found for the real samples. 

Since every burst is counted only once, the maximum 
possible number of coincidences is equal to the number of 
GRBs used. This method is justified since any observed 
burst can have only one counterpart. With this method of 
counting only once, the obtained number of coincidences 
depends both on the number of objects in the compari- 
son sample as well as on their distribution on the sky, in 
particular their isotropy. In Fig. 1 the value < Cb > is 
plotted against the size of the various comparison sam- 
ples of QSOs/AGNs we studied (see Sect. 3 for details). 
Clearly, Cb is not a linear function of the size of the com- 
parison sample. A linear correlation holds for most of the 
samples with less than approximately 2000 members. For 
larger samples the < Cb > values obtained are signifi- 
cantly lower than expected from a linear extrapolation of 
the smaller samples. However, with the 1000 simulations 
of the GRB positions we establish the distribution of the 
coincidences for a given sample. This enables us to ob- 
tain the corresponding probability for a real correlation of 
GRBs with the respective comparison sample. 
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2.3. Shifted GRB sample 

The small-scale structure in the GRB locations is more 
problematic than the large-scale structure. Although no 
evidence for an autocorrelation in the 3B catalogue has 
been reported, a possible intrinsic correlation of GRB po- 
sitions can not be excluded at present, given their large 
positional uncertainties. Rather than testing the 3B cat- 
alogue itself for the presence of small-scale structure, we 
will test our results using a method intrinsically robust 
against a possible hidden small-scale structure in the GRB 
distribution. 

We applied the same method as in the previous section 
to calculate the number of coincidences. However, a differ- 
ent method to construct the simulated random positions 
was introduced to determine the number of chance coin- 
cidences. One thousand control samples were generated 
by changing the sign of the declination of the measured 
burst positions and shifting the right ascension 1000 times 
with a random number in the range between 5° and 255°. 
The random samples constructed in such a way conserve 
the small-scale structure of the real GRB distribution and 
allow to test its influence on the correlations. 



2.4- Singular Coincidences 

By using the shifted GRB sample we are able to test our 
results against a possible small-scale structure of the GRB 
position distribution. Because the QSO/AGN positions 
provided in the Veron-Cetty & Veron (1996) compilation 
are not distributed homogeneously over the sky, we need 
a method which is robust against "clustering" of objects 
in the comparison samples. This is especially important 
as most of the "clustering" is due to the accumulation of 
objects in several deep small-area searches for AGN and 
QSOs. To this end we counted only bursts with a single 
QSO or AGN in their 1 a error circle. With this method 
we exclude every burst with two or more QSO / AGN coun- 
terparts, i.e. too many to allow a direct comparison with 
the large error boxes of GRBs. The expected background 
values and the corresponding expected distribution for an 
uncorrelated burst sample were obtained again with ran- 
domly simulated GRB samples. 

The statistical errors of the GRB positions range from 
0?13 to 34?15. Including the systematic error this results 
in a range from 1?61 to 34? 18 for the total error. The mean 
total error for all bursts is 5?2±4?1, and the mean area of 
the error circles is 138 deg 2 per burst. Therefore, the sum 
of all 68% confidence error boxes of 3B bursts corresponds 
to more than three times the whole sky. In order to avoid 
too many spurious coincidences, we restricted our study to 
the bursts with small errors. In Table 1 we list the charac- 
teristics of the different GRB samples used here. As men- 
tioned above, both the statistical and total errors depend 
on the Poisson uncertainty in the BATSE measurements. 
Therefore, bursts with higher flux generally show smaller 



positional errors. In Fig. 2 the total error radii for 835 
bursts are shown versus the burst flux in the 20-50 kcV 
energy band. Thus, restricting the bursts to those with 
a smaller error circle also restricts them to a higher flux 
range. If we search for positional coincidences of GRBs 
with AGNs and quasars and thus assume that GRBs orig- 
inate at cosmological distances, then the restricted flux 
range may correspond to a restricted redshift range. 

Table 1. Characteristics of the GRB Samples Used 



Maximal Error 


Number 


total 


statistical 


of bursts 


Rt < 1°7 


r 3t < 0?574 


80 


Rt < 1?8 


r st < 0?825 


134 


Rt < 2?0 


r st < 1?200 


206 



2.5. The AGN /Quasar Samples Studied 



We used the latest (7 th ) version of the "Catalogue of 
Quasars and Active Nuclei" (Veron-Cetty & Veron 1996). 
This catalogue has been compiled from the available lit- 
erature and provides the position of 11,662 sources drawn 
from 1662 references. Due to several deep small-area 
searches for AGN and QSOs the objects in this compila- 
tion are not distributed homogeneously over the sky, but 
appear "clustered" in these areas. Therefore, any test for a 
positional correlation between GRBs and quasars or AGN 
must be robust against observational bias inherent in the 
analyzed sample. 
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Fig. 2. The total position error is shown over the measured 
flux in the 20-50 keV band for 835 bursts. The total position 
error is always greater than the systematic error of 1?6. 
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Table 2. Comparison samples of QSO/AGN and positional coincidences with three GRB samples of different positional errors. 
< Cb > is the expected number of chance coincidences derived from 1000 randomly simulated GRB samples. Cm gives the 
number of coincidences with the 3B catalog bursts and P(%) is the fractional number of simulations revealing less coincidences 
than Cm- 



Main Subclass N< a > D {b) R t =l°7 R t =l°8 R t =2°0 

Class or (°) 

selection < C B > C 1 o 7 /P(%) < C B > C 1 o g /P(%) < C B > C f JP{%) 



QSO 


all 


8609 


0.6 ±0.9 


31.3 ±4.2 


QSO 


radio- quiet 


7146 


0.4 ± 0.9 


20.3 ± 3.9 


QSO 


radio-loud 


1377 


2.3 ± 1.6 


18.2 ± 3.6 


QSO 


highly polarized 


72 


10. ±7.4 


1.2 ± 1.1 


BL 


all 


220 


5.9 ±4.3 


3.5 ± 1.9 


BL 


confirmed 


93 


8.6 ±5.2 


1.6 ± 1.2 


BL 


highly polarized 


76 


11. ±6.8 


1.3 ± 1.1 


BL 


radio selected 


119 


8.4 ±5.5 


1.9 ± 1.4 


BL 


X-ray selected 


82 


8.7 ±5.0 


1.4 ± 1.2 


AGN 


AGN 


1553 


2.1 ± 1.7 


19.0 ±4.0 


AGN 


Seyfert 1 


888 


2.7 ±2.2 


12.1 ± 3.4 


AGN 


Seyfert 2 


496 


3.9 ±2.6 


7.3 ±2.6 


AGN 


Prob. Seyfert 


97 


8.0 ±6.1 


1.4 ± 1.2 


AGN 


LINER 


71 


10. ± 7.8 


1.2 ± 1.1 


AGN 


radio-quiet 


1346 


2.2 ± 1.8 


16.7 ±3.8 


AGN 


radio-loud 


166 


7.4 ±4.2 


2.8 ± 1.6 




AG (C) 


1165 


1.0 ±2.1 


5.8 ±2.3 




Nucl. HIT 


116 


7.7 ±5.9 


1.8 ± 1.3 



36 


/ 


83. 


.4 


53 


6±5 


.7 


63 


/ 94. 


2 


85 


.3 ±7 


.0 


98 


/ 


95.2 


28 


/ 


96. 


.4 


34. 


7±5 


.1 


44 


/ 94. 


6 


55 


.3 ±6 


.4 


66 


/ 


94.0 


15 


/ 


16 


3 


31 


5±4 


.9 


35 


/ 73. 


1 


51 


,1 ±6 


.2 


58 


/ 


84.8 


1 


/ 


31. 


.5 


2. 


0±1 


.4 


3 


/ 67, 


8 


3 


,4± 1 


,7 


5 


/ 


74.8 


2 


/ 


13 


,5 


6 


2±2 


.5 


3 


/ 4. 


7 


10 


.3 ±3 


.2 


5 


/ 


2.0 


2 


/ 


52. 


9 


2. 


8± 1 


.6 


2 


/21. 


6 


4 


,5 ±2 


.1 


3 


/ 


16.4 


1 


/ 


27. 


.6 


2 


2± 1 


.5 


2 


/ 35. 


.1 


3 


,7± 1 


.9 


3 


/ 


28.7 





/ 








3 


4± 1 


.9 


1 


/ 3. 


5 


5. 


.7 ±2 


.4 


2 


/ 


2.7 


2 


/ 


59 


.5 


2. 


4± 1 


.5 


2 


/31. 


2 


4 


,0±2 


.0 


3 


/ 


24.0 


20 


/ 


55. 


.1 


32. 


8±5 


.2 


34 


/ 55, 


5 


53 


,2 ±6 


.4 


55 


/ 


57.4 


13 


/ 


55. 


.7 


20. 


.8 ±4 


.3 


21 


/ 47, 


.7 


34 


,0±5 


.4 


36 


/ 


64.3 


6 


/ 


26. 


.2 


12. 


.7 ±3 


.5 


12 


/ 38, 


2 


21 


,0 ±4 


.4 


20 


/ 


38.2 


3 


/ 


82. 


.5 


2. 


6± 1 


.6 


3 


/ 54, 


6 


4 


,3 ±2 


,1 


3 


/ 


19.5 





/ 





.0 


2 


.1 ± 1 


.4 


1 


/14. 


8 


3 


,4± 1 


.9 


2 


/ 


14.1 


18 


/ 


58. 


6 


28. 


.8 ±4 


.9 


31 


/ 62, 


.8 


46 


.8 ±6 


,1 


50 


/ 


66.9 


1 


/ 


6 


.7 


4. 


8±2 


.1 


1 


/ o. 


.4 


7. 


.9 ±2 


.8 


2 


/ 


0.3 


9 


/ 


87. 


.5 


10. 


.1 ±3 


.0 


18 


/ 98, 


.7 


16 


.3 ±3 


.9 


25 


/ 


98.1 


1 


/ 


15. 


,7 


3 


2± 1 


.7 


1 


/ 3, 


6 


5 


,3 ±2 


,2 


2 


/ 


2.6 



( a > Number of sources. 

^ Actual mean distance to nearest neighbour of the sample and its standard deviation. Since the sources in these samples are 
distributed rather irregularly, the standard deviation occasionally exceeds the mean. 

Active Galaxies are galaxies with unusually high X-ray or radio (few) luminosity. These do not overlap with any other 
subsample in this table (see text). 



The "Catalogue of Quasars and Active Nuclei" (Veron- 
Cetty & Veron 1996) is divided in three parts, — quasars, 
BL Lac objects and AGN. In addition to these three sub- 
samples we shall consider certain subsamples, e.g. sources 
selected by spectral type. Because synchrotron emission 
is a favoured mechanism to explain GRBs, we used crite- 
ria which are thought to favour the presence of jets. The 
resulting samples are listed in Table 2 together with the 
sample size and the mean distance to the nearest neigh- 
bour. The selection criteria for the subsamples chosen are 
described in the following. 

2.5.1. Quasars 

We adopt the definition by Veron-Cetty & Veron (1996) 
of a quasar as a star-like object or an object with a 
star-like nucleus, brighter than absolute visual magnitude 
—23. In addition to the total quasar sample we consider 
three subsamples. The first is given by objects charac- 
terized by a high optical polarization, P > 3%, which is 
mostly interpreted as evidence for synchrotron radiation 
connected with a Doppler-boosted nuclear jet (e.g. Wills et 
al. 1992). The other two classes are defined by their radio- 
loudness. Although the radio energy band is far away from 



the gamma-ray band discussed here, these subsamples are 
justified because of the strong correlation found between 
the source properties in the radio and X-ray region. Al- 
ready Einstein data showed a correlation between radio- 
loudness and both X-ray luminosity and X-ray spectral 
index (Wilkes & Elvis 1987). This result was confirmed 
by ROSAT measurements based on larger samples (see 
e.g. Schartel et al. 1996). In particular, the different spec- 
tral properties may indicate that the X-ray spectra are 
dominated by different physical mechanisms which may 
result in a different behavior at higher energies. We define 
the radio-loudness based on the radio-to-optical spectral 
index (Stocke et al. 1985; Kellermann et al. 1989) 

1 , / 55GHz \ 

aRO - K38 l0S I ^ T ) 

\ 2500A / 

= 5~b [ log ^ 5GHz ~ 23 -° + 19 - 76 + 0AV 

-0.3 log (z + 1)] 

where the spectral fluxes are given in units of Jy and V 
is the object's apparent visual magnitude. Note that this 
definition is independent of optical luminosity. 
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The radio-to-optical index has been K-corrected as- 
suming a spectral index of an — —0.7 in the radio range 
and ao = — 1.0 in the optical regime. The radio fluxes 
and the visible magnitudes were taken from Veron-Cetty 
& Veron (1996). We define radio-loud quasars as those 
with oibo > 0.3. A source is defined to be radio-quiet if 
either a^o < 0.3 or no radio flux is given in the catalogue. 

2.5.2. BL Lac objects 

The Veron-Cetty & Veron (1996) compilation provides 
225 confirmed, probable or possible BL Lac objects. In ad- 
dition to the total sample and the sample restricted to the 
optically confirmed BL Lacs we consider three further sub- 
samples. The first one contains all known BL Lac objects 
which show an optical continuum polarization P >3% and 
the second one (called "radio-selected" ) is given by all ob- 
jects with a flux density > O.Uy at either A = 6 or 11 cm. 
Finally, we select all sources not fulfilling one of the above- 
mentioned criteria. This latter sample is dominated by BL 
Lac objects which are strong X-ray emitters. For 66 of the 
82 objects in this group the mere designation already in- 
dicates that they were discovered during X-ray missions. 

2.5.3. AGN 

The Veron-Cetty k Veron (1996) compilation lists 2833 
AGN. By definition of the authors, all of them are 
fainter than absolute magnitude -23. This sample in- 
cludes sources with spectra classified as Seyfert 1, Seyfcrt 
2 or LINERs following the definition of Heckman (1980). 
Further 97 sources are classified as probable or possible 
Seyfert galaxies. 

The list of AGN includes 116 galaxies with a nuclear 
HII region. These are now interpreted as a burst of star 
formation, but were called Seyferts in the past and later 
were reclassified. Therefore, we treat them separate from 
the Seyferts here. We collect all objects without any mor- 
phological type assigned to them by Veron-Cetty and 
Veron (1996) into a class we term "Active Galaxies" (here- 
inafter AG). These typically consist of galaxies with nu- 
clear radio or X-ray sources. Note that we use the term 
"AGN" only for the combined sample of Seyfert galax- 
ies and LINERs. We also divide the AGN into radio-loud 
and radio-quiet samples, defined by their radio-to-optical 
index aoR, in the same way as above for the quasars. 

3. Search for Coincidences 

3.1. Positional Coincidence Analysis 

In order to test our results if only for internal consistency, 
we compare the subsamples of QSOs and AGN defined 
above with three samples of GRBs containing bursts up 
to a given positional error radius (see Table 1). The results 
are listed in Table 2. For every comparison we give the ex- 
pected chance coincidences < Cb > derived from 1000 



randomly simulated GRB samples. Further we character- 
ize the distribution of C by its standard deviation. Since 
the distribution of coincidences is not Gaussian, we also 
give the number of coincidences obtained with the real 
GRB samples, C 1 o ? , C 1 o & , and C 2 o Q , together with the 
probability P to measure a lower value than the one ob- 
tained. Histograms of the distribution of C, i.e. the num- 
ber of coincidences for the artificial sets of GRBs (null hy- 
pothesis), are shown in Fig. 3 for all analyzed subclasses. 
The expected value of coincidences (or the mean of the 
distribution), < Cb >, is shown as dotted line and the 
obtained value of real coincidences, C 1 o ? is indicated by a 
dashed line. 

Three classes, the "active galaxies", all quasars and 
the radio-quiet quasars show coincidences in excess of the 
expected value at the 90%-level or more. Since the radio- 
loud quasars alone do not show any excess coincidences, 
we conclude that only radio-quiet quasars may be related 
to GRBs. All other classes are compatible with the null 
hypothesis of no relation with GRBs. In particular, the 
highly-polarized and the radio-loud subsamples show no 
significant probability of being associated with GRBs. 

3.2. Radio-quiet Quasars and Active Galaxies 

Only two classes, active galaxies and radio-quiet quasars 
show coincidences in excess of random expectation above 
the 95%-level. The sample sizes range from 71 to 7146 
members and the mean distance to the nearest neigh- 
bour within a given sample varies from 0?4 to 11° (see 
Tablc2). In this section we will test whether any excess 
correlation may be influenced by the size of the compari- 
son sample and/or the distance to the nearest neighbour. 
The expected number of chance coincidences was obtained 
by comparing the positions of the sources with randomly 
simulated positions of the GRBs (as explained in Sect. 
2). Therefore, our results are robust against large-scale 
structure in the burst locations. A possible influence of 
any small-scale structure in the burst locations remains 
to be checked. We will also study the effect of the non- 
homogeneous distribution of objects in the comparison 
samples. In summary, we have to check for any hidden, 
yet unreported characteristics in the sky distribution of 
GRBs and/or in the sky distribution of AGN and QSOs 
which may be responsible for the high levels of correlation 
with GRBs. 

As mentioned above the position errors of GRBs are a 
function of the measured 7-ray flux above 25 keV. A re- 
striction to bursts with the smallest errors leaves us with 
the brightest bursts and thus more coincidences for nearby 
GRB sources are expected. Therefore, sample selection by 
redshift offers an additional possibility to test our result, 
as we would expect more coincidences for the low-redshift 
sources. The main advantage of redshift-selected subsam- 
ples is given through their intrinsic parameters (for ex- 
ample the mean distance to the next source) which are 
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Fig. 3. Distribution of the position coincidences for 1000 random simulations of the 80 burst positions with R <1?7. The dotted 
line shows the expected value for the null-hypothesis and the broken line the number of real coincidences between the gamma-ray 
burst positions and positions of quasars and AGN defined by different samples: The distribution is given for: quasars in the first 
panel: (a) all, (b) radio-quiet, (c) radio-loud, and (d) highly polarized quasars, for BL Lac objects in the second panel: (e) all, 
(f) confirmed, (g) highly polarized, (h) radio-selected, and (i) X-ray selected, for AGN in the third and fourth panel: (j) all, (k) 
Seyfert 1, (1) Seyfert 2, (m) probable Seyferts and (n) Liners, (o) radio-quiet, (p) radio-loud. Finally, the distribution is shown 
for active galaxies (s) and Nuclear H II regions (r). A detailed description of the selection criteria is given in the text. 
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Table 3. Positional coincidences of redshift subsamples with randomly simulated GRB samples. 



Redshift 


Number 


Distance to 


Rt- 


=1?7 


Rt 


=1?8 








of 


next Source 


<C B > 


C 1?7 /P(%) 


<C B > 




/P{%) 




_ 

z 


sources 


D(°) 
















Radio-quiet Quasars 








0.00—0.58 


0.35 


883 


2.1 ± 2.1 


9.7 ±3.0 


17 / 98.1 


16.8 ± 3.9 


28 


/ 99.5 


0.58 — 1.03 


0.81 


889 


1.3 ± 1.7 


6.2 ±2.4 


13 / 99.1 


10.8 ± 3.2 


22 


/ 99.8 


1.03—1.41 


1.24 


903 


1.1 ± 1.9 


5.5 ±2.3 


9 / 90.0 


9.6 ±3.1 


12 


/ 73.4 


1.41—1.73 


1.57 


895 


1.1 ± 1.8 


5.4 ±2.2 


8 / 82.7 


9.3 ±2.9 


11 


/ 66.4 


1.73—2.00 


1.89 


871 


1.2 ±2.2 


5.8 ±2.2 


7 / 66.8 


9.9 ±2.9 


12 


/ 73.0 


2.00—2.21 


2.11 


911 


1.3 ±2.2 


6.0 ±2.3 


8 / 75.9 


10.5 ±3.0 


11 


/ 50.9 


2.21—2.50 


2.33 


887 


1.2 ± 1.8 


6.1 ± 2.3 


6 / 41.8 


10.5 ±3.1 


12 


/ 63.0 


z > 2.50 


2.98 


881 


1.4 ±2.0 


6.8 ±2.6 


8 / 63.2 


11.9 ±3.3 


16 


/ 85.5 








Active Galaxies 










z < 0.17 


0.07 


590 


1.4 ±3.4 


3.5 ± 1.8 


6 / 87.6 


6.1 ±2.4 


12 


/ 97.4 


z > 0.17 


0.40 


575 


1.5 ±2.7 


4.2 ±2.0 


7 / 87.0 


7.2 ±2.6 


14 


/ 98.9 



different from the total sample. Therefore, they allow to 
test our results against these parameters. Table 3 lists the 
redshift bins, their mean redshift, the number of sources in 
each bin, and the mean distance to the nearest neighbour. 



3.2.1. Distance to the nearest neighbours 

Table 2 shows that comparison samples with the high- 
est correlation with GRBs (with the exception of the to- 
tal quasar sample) have the smallest mean distance to 
the nearest neighbour. Active galaxies have a mean dis- 
tance of 1?0±2?1 to their nearest neighbour and radio- 
quiet quasars have only 0?4±0?9. The subsamples defined 
by redshift allow to increase significantly the distance to 
the nearest neighbour, as is obvious from Table 3. In par- 
ticular, in the lowest redshift sample of radio-quiet quasars 
a mean distance of 2?1±2?1 is achieved. Samples with 
comparable values of nearest neighbour distance, e.g. all 
AGN (2?1±1?7), radio-quiet AGN (2?2±1?8) and radio- 
loud quasars (2?3±1?6) do not show coincidences with 
GRBs in excess of random expectation. We correlated 
the redshift-selected samples of QSO/AGN with randomly 
simulated versions (see Sect. 2) of the two GRB samples 
with R t < 1?7 and R t < 1?8. The results are given in 
Table 3. For the two subsamples of radio-quiet quasars 
with lowest redshift, probabilities of > 98% are found for 
a correlation with both GRB samples. Such a positive re- 
sult is not obtained for the "active galaxies" , which show 
probabilities as low as 87% for the GRB sample with the 
smallest error radii. 

The z < 0.58 subsample for the radio-quiet quasars 
shows that our result is robust against an increase of the 
mean distance to the nearest neighbour. On the other 
hand, the highcr-redshift subsamples demonstrate that 
even samples with a small distance to the nearest neigh- 
bour may well yield negative results. 



3.2.2. Small-scale structure in the GRB distribution 

In order to test the positive result obtained for radio-quiet 
quasars and to investigate in more detail the results ob- 
tained for active galaxies for a possible bias by the small- 
scale structure of the GRB distribution, we tested our re- 
sults using shifted GRB samples (cf. Sect. 2.3). This covers 
two different effects: First, local differences in the BATSE 
exposure time of certain sky regions causing fluctuations 
of the local GRB density are smeared out. Second, any 
possible internal structure of the GRB distribution (e.g. 
similar to the Great Wall structure of galaxies) is tested. 

The results are presented in the usual way in Table 4. 
They strengthen the results obtained in the previous sec- 
tion. In the first two redshift subsamples (z < 0.58 and 
0.58 < z < 1.03) the radio-quiet quasars reach proba- 
bility levels of >99% for correlation with GRBs for both 
burst samples. In contrast, the two subsamples of the ac- 
tive galaxies do not exceed a probability of 94% for corre- 
lation with GRBs. 



Table 4. Redshift Subsamples / Shifted comparison samples 



Redshift 


Rt 


=1?7 




Rt 


=1?8 






Range 


<C B > 


C 1?7 /P(%) 


<C B > 


C l°8 


/P(%) 




Radio-quiet Quasars 










0.00—0.58 


9.6 ±2.8 


17 / 99.8 


16.6 ± 


3.9 


28 


/ 


99.7 


0.58—1.03 


6.0 ±2.0 


13 / 99.9 


11.0 ± 


2.9 


22 


/ 


100. 


1.03—1.41 


5.2 ± 2.0 


9 / 94.2 


9.6 ± 


2.9 


12 


/ 


74.9 


1.41—1.73 


4.8 ± 2.1 


8 / 89.6 


8.9 ± 


3.0 


11 


/ 


69.3 


1.73—2.00 


5.2 ± 2.1 


7 / 70.9 


9.5 ± 


3.1 


12 


/ 


74.7 


2.00—2.21 


5.5 ±2.2 


8 / 81.2 


9.9 ± 


3.3 


11 


/ 


59.7 


2.21—2.50 


5.2 ± 2.3 


6 / 57.7 


10.1 ± 


3.2 


12 


/ 


67.9 


z > 2.50 


6.1 ±2.5 


8 / 73.6 


11.6 ± 


3.1 


16 


/ 


87.9 






Active Galaxies 










z < 0.17 


3.6 ± 1.9 


6 / 83.8 


7.6 ± 


2.3 


12 


/ 


93.4 


z > 0.17 


4.2 ± 2.0 


7 / 88.8 


8.4 ± 


2.7 


14 


/ 


93.7 
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Table 5. Redshift Subsamplos / Singular Coincidences 



Kedsnut 


Rt 


=1?7 


Rt- 


1 ° O 

=1. 8 


Range 


<C B > 


C 1?7 /P(%) 


<C B > 


C 1?8 /P(%) 




Radio-quiet Quasars 




0.00—0.58 


6.9 ±2.6 


14 / 99.0 


11.8 ±3.4 


23 / 99.6 


0.58—1.03 


3.7 ± 1.9 


8 / 95.9 


6.4 ±2.5 


15 / 100. 


1.03—1.41 


2.9 ± 1.7 


3 / 43.8 


5.1 ± 2.2 


4 / 24.9 


1.41—1.73 


2.9 ± 1.7 


1 / 6.0 


5.0 ±2.2 


3 / 11.1 


1.73—2.00 


3.1 ± 1.7 


4 / 63.2 


5.4 ±2.3 


7 / 72.4 


2.00—2.21 


3.3 ± 1.7 


3 / 32.8 


5.9 ±2.3 


5 / 28.4 


2.21—2.50 


3.5 ± 1.7 


2 / 12.4 


6.0 ±2.4 


5 / 27.0 


> 2.50 


4.0 ±2.0 


5 / 62.8 


7.0 ±2.6 


8 / 61.1 






Active Galaxies 




z < 0.17 


1.7 ± 1.3 


2 / 46.8 


3.0 ± 1.7 


4 / 65.8 


z > 0.17 


2.5 ± 1.5 


3 / 54.2 


4.3 ±2.0 


6 / 74.7 



3.2.3. Non-Homogeneity of the Comparison Samples 

In Sect. 3.2.1 we were able to rule out a bias caused by 
the mean distance to the nearest neighbours for the radio- 
quiet quasars whereas our results for the active galax- 
ies were not convincing. We are aware that the non- 
homogeneity of the comparison sample is not sufficiently 
characterized by the mean distance to the nearest neigh- 
bours. To test for a possible bias of this effect, we count 
only bursts with only a single QSO or AGN in its 1 a 
error circle (see singular coincidences in Sect. 2.4). The 
expected chance coincidences and their distribution were 
obtained again with simulated GRB positions. 

Table 5 gives the results and confirms a high probabil- 
ity for correlation with GRBs in excess of random expec- 
tation for the two lowest redshift quasar samples, reaching 
a probability above the 99% level. In contrast to this, we 
do not find any support for a correlation with GRBs for 
the "active galaxies" . 

4. Results 

We first summarize all classes for which no coincidences in 
excess of random expectation were found. In the second 
subsection we will strengthen the case for a correlation 
between GRBs and radio-quiet quasars. 

4-1- Upper Limits 

Quasars and AGN are often discussed in the frame of 
the standard model assuming a massive compact object 
in the center, most probably a super-massive black hole. 
Many spectral features, like optical polarization or radio- 
loudness, are interpreted as evidence of a nuclear jet. Most 
of the samples and subsamples of QSO and AGN analyzed 
here show no evidence for a correlation with GRBs. How- 
ever, the measured values can be used to estimate upper 
limits for the emission of GRBs. To determine the 3 a up- 
per limit we use always our analysis ± the R t < 1?8 burst 



Table 6. Upper Limits 



Main 


O LI LtUldSS 


1\I 1 1 m ripv 
IN UlilUtl 






oiass 


or 


OI 






selection 


sources 






QSO 


radio- loud 


1377 


18.2 


1.50 


QSO 


highly polarized 


72 


5.2 


8.17 


BL 


all 


220 


4.3 


2.21 


BL 


confirmed 


93 


4.0 


4.87 


BL 


highly polarized 


76 


4.3 


6.40 


BL 


radio-selected 


119 


3.3 


3.14 


BL 


X-ray selected 


82 


4.1 


5.66 


AGN 


AGN 


1553 


16.8 


1.22 


AGN 


Seyfert 1 


888 


13.1 


1.67 


AGN 


Seyfert 2 


496 


9.8 


2.24 


AGN 


Probably Seyfert 


97 


5.2 


6.07 


AGN 


LINER 


71 


3.1 


4.94 


AGN 


radio-quiet 


1346 


16.9 


1.42 


AGN 


radio-loud 


166 


2.5 


1.70 




Active galaxies (AG) 


1165 


16.9 


1.64 




Nuclear H II region 


116 


2.9 


2.83 



: per source in 10 2 yr 



sample which reaches the highest significance. The 3 a up- 
per limit is given by: 

Us a = CjOg ± SACjOg— < CjOg > 

where C 1 o s is the number of coincidences with the real 
GRBs, and < C 1 o s > is the number of coincidences ex- 
pected for a simulated burst sample. ACjOg is approxi- 
mated by the standard deviation of < C 1 o g > as provided 
in Table 2, together with the other two values. The re- 
sulting upper limits are listed in Table 6. As a quantity 
independent of the amount of members in the different 
subclasses we define the emission rate, Ru 3a ■ 

R U3a = (1/0.38) (1/0.68) (1/3.42) U^N' 1 per year 

= 1.132 U^N^ 1 per year . (2) 

This rate is given by the derived upper limits divided by 
the size N of the comparison sample, and by the observa- 
tion time of BATSE (3.42 years). Further it is corrected 
for the 68% confidence interval we used, and the proba- 
bility of detecting a burst (38%, cf. Meegan et al. 1996). 
All numbers refer to the 134 GRBs in the 3B catalogue 
having R t < 1?8. 

Our initial positive result for a correlation of GRBs 
with "active galaxies" (AG) was not confirmed in the pre- 
vious chapter. Thus we give an upper limit for AGs in 
Table 6. The positive result obtained initially might be 
due to the non-homogeneity of the sky distribution of AG. 
However, an additional bias due to a small-scale structure 
of the GRB distribution can not be excluded at this point. 
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Table 7. Magnitude subsamples of radio-quiet quasars 



Magnitude Number R t =l°7 R t =l°8 

of 



Range 


M abs 


sources 


<C B > 


C 1?7 /P(%) 


<C B > 




-23.0 24.1 


-23.6 ±0.3 


836 


7.3 ±2.6 


9 / 67.4 


12.6 ±3.3 


18 / 92.4 


-24.2 25.0 


-24.6 ± 0.3 


847 


5.9 ±2.3 


13 / 99.4 


10.2 ±3.0 


18 / 99.2 


-25.1 25.6 


-25.4 ±0.2 


784 


5.2 ± 2.2 


9 / 91.9 


8.9 ±2.9 


12 / 81.9 


-25.7 26.1 


-25.9 ±0.1 


820 


4.8 ± 2.2 


8 / 88.0 


8.3 ±2.8 


15 / 98.2 


-26.2 26.6 


-26.4 ±0.1 


910 


5.7 ±2.2 


9 / 90.2 


9.7 ±3.0 


17 / 97.9 


-26.7 27.1 


-26.9 ±0.1 


875 


6.1 ±2.3 


11 / 96.8 


10.6 ±3.1 


16 / 93.3 


-27.2 27.9 


-27.5 ±0.2 


993 


7.3 ±2.6 


10 / 81.4 


12.6 ±3.3 


15 / 71.9 


-28.0 31.9 


-28.8 ±0.7 


871 


8.9 ±3.0 


7 / 21.1 


15.6 ±3.7 


14 / 29.9 



4-2. Radio- quiet Quasars 

Because of the inverse relation between burst flux and po- 
sitional error radius we expected to strengthen the result 
obtained in our first search for correlation for the low- 
rcdshift quasars. In Sect. 3.2. we were able to show that 
the result is influenced neither by the distance to the near- 
est neighbour, nor the small-scale structure of the burst 
locations, nor by the non-homogeneity of their sky dis- 
tribution. In this section we will try to strengthen our 
positive correlation for radio-quiet quasars. 

Seyfert galaxies are closely related to radio-quiet 
quasars. Because many radio-quiet quasars show a Seyfert- 
like spectrum, it appears that the main difference between 
these two classes is the absolute magnitude. Given the 
negative result obtained for the Seyfert galaxies, it seems 
reasonable to presume that the absolute magnitude in- 
fluences the emission rate of GRBs. Therefore, we also 
defined subsamples according to absolute magnitude as 
listed in the compilation of Veron-Cetty & Veron (1996). 
Again, the values obtained were compared with randomly 
simulated GRB samples and the results are provided in 
Table 7. The probability levels do not support coincidences 
in excess of random expectation for the fainter magnitudes 
(-23.0 > M abs > -24.1) at more than the 95%-level. 
However, excess coincidences are indeed found for quasars 
brighter than M a \, s = —24.2. It appears that only intrin- 
sically bright quasars show a detectable emission rate of 
GRBs. This interpretation is in agreement with the nega- 
tive result found for the Seyfert galaxies. 

We conclude that there is an increased probability 
for the nearby and for the absolutely brighter radio-quiet 
quasars to correlate with GRBs. It is therefore tempting to 
construct an "optimized quasar sample" with even higher 
probability for emitting GRBs. We defined such a sample 
by z < 1.0 and M a b s < —24.2. It contains 967 radio- 
quiet quasars and was tested with the three methods de- 
scribed in Sect. 2, for the two GRBs samples with the 
smaller error radii. Results are provided in Table 8. The 



Table 8. Quasar sample selected by magnitude and redshift 



Method" Rt=\°7 R t =l°8 

<C B > C 1 o 7 /P{%) <C b > C 1?g /P(%) 

Radio-quiet Quasars with : 
Magnitude < -24.2 and z < 1.0 

1 8.4 ± 2.8 20 / 100.0 14.5 ± 3.7 27 / 99.7 

2 7.9 ± 2.3 20 / 100.0 14.0 ± 3.5 27 / 100.0 

3 5.4 ± 2.3 17 / 100.0 9.2 ± 3.1 20 / 99.8 

a Numbers denote the different samples: 1 = simulated control 
sample; 2 = shifted control sample; 3 = singular coincidences 



corresponding distribution of the positional coincidences 
for the null-hypothesis of the "optimized quasar sample" 
are plotted in Fig. 4. The value expected at random is 
shown as a dotted line and the value obtained for the 
real GRBs is indicated by a dashed line. For all six com- 
parisons the probability for a real correlation with GRBs 
exceeds 99.7%. Moreover, in four of the six comparisons 
the 1000 simulations of GRBs positions did not show a 
single case with a number of coincidences as high as for 
the real GRBs. 

Frequently, two catalogues of objects detected in dif- 
ferent wavebands are cross-identified from a study of the 
radial histogram of the distances between suspected coun- 
terparts. We did not use this method until this point be- 
cause many bursts have more than one QSO or AGN in 
their error circle, which would significantly increase the 
chance coincidences at large distances between GRBs and 
objects of the comparison sample. However, the "opti- 
mized quasar sample" (z < 1.0 and M a b s < —24.2) allows 
to illustrate our result with the radial histogram of the 
distances between suspected counterparts. The differen- 
tial distribution of coincidences is given by the number of 
quasars, dN in an annulus between a radius r and r ± dr. 
The distribution of the real coincidences is a function of 
the positional errors which is not very well determined. 
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Fig. 4. Distribution of the position coincidences of the null-hypothesis for the optimized quasar sample given by z < 1.0 and 
M a i,s < —24.2. The expected value of the background position coincidences is indicated with a dotted line and the obtained 
value of the real coincidences is shown with a dashed line. The upper panel shows the comparison for 80 gamma-ray burst 
positions defined by R<1?7 and the lower for 132 burst positions defined by R<1?8, respectively. In a and d the 80 (132) 
positions were simulated randomly 1000 times. The distributions b and e were obtained by shifting the original burst positions 
1000 times. And finally, in c and f only bursts are counted which show not more than one coincidence. 



However, the real coincidences must be located within a 
distance less than 3 times the positional error radius, R t , 
assuming a two-dimensional Gaussian distribution. The 
number of randomly expected coincidences is proportional 
to the cumulative area of all annuli. Thus, the number of 
sources in a given annulus as a function of the radius is 
given by 

dN(r) = A{r, AD)2irrdr + B2irrdr 

where A(r) is the probability for real coincidences per unit 
area, and B is the probability for chance coincidences per 
unit area. We divide by r and set 

C = A/(2n) and D = B/(2n). 



Thus, we may rewrite: 

= C(r, AD)dr + Ddr 

r 

To allow for the very different sizes of the GRB po- 
sitional errors we scale the distance between GRB and 
radio-quiet quasars in units of the positional error a, i.e. 
we choose dr = a to obtain the normalized differential dis- 
tribution, which should be constant for r > 3cr. In Fig. 5 
the normalized differential distribution for the optimized 
quasar sample is given. In the normalized differential dis- 
tribution every coincidence between a gamma-ray burst 
and a quasar is counted. This results in a doubling of the 
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Fig. 5. The area-normalized distribution of pairs of gamma-ray 
bursts and radio-quiet quasars as function of their mutual dis- 
tance, scaled in units of the position uncertainty a of the 
bursts. The left panel shows only the N=80 GRBs with po- 
sition uncertainty of Rt < 1?7. The right panel is based on 136 
GRBs with R t < 1?8. In both figures the sample of N = 967 
radio-quiet quasars is defined by z < 1.0 and M a bs < —24.2. 

background value for the first bin as opposed to counting 
every GRB only once. Nevertheless, the coincidences we 
find in excess of random expectation between 2 a and 5 a 
more than compensate for the poor statistical significance 
of the first bin (cf. Fig. 5). 

By substituting {7 3(T with (C^Og— < C 1 o g >) in equa- 
tion 2 we obtain the emission rate for the detected radio- 
quiet quasars: 

Remission = 1-132 (C^g- < C 1?8 >)N^ in y" 1 

Using the results obtained for the " optimized quasar sam- 
ple" given in the first line of table 8 we determine an 
emission rate of burst of 0.0146±0.004 y _1 per radio-quiet 
quasar. 

5. Discussion 

5.1. Correlation with radio-quiet quasars 

We have found evidence that intrinsically bright radio- 
quiet quasars with z < 1.0 coincide with GRBs from the 
3B catalogue at a confidence level of more than 99%. 

The correlation with low-redshift objects is not surpris- 
ing. For statistical reasons our study had to be restricted 
to GRBs with the smallest positional errors, which tend to 
be the brightest GRBs. These GRBs are preferably nearby 
since the luminosity distribution of GRBs cannot be very 
flat (Hakkila et al. 1995), and thus one introduces a strong 
selection effect for low-redshift quasars. Based on the dif- 
ferent methods which we applied to test our results, we 
are confident to have excluded any bias from an inhomoge- 
neous sky distribution of quasars. In addition, we showed 
our result to be robust against a possible small-scale struc- 
ture in the sky distribution of 7-ray bursts. 

The restriction to intrinsically bright quasars with 
M a bs < —24.2 seems to be a real physical selection indicat- 



ing that the GRB emission rate is at least lower for intrin- 
sically faint quasars. This interpretation is supported by 
the negative result obtained for the Seyfert galaxies which 
are intrinsically fainter than radio-quiet quasars but oth- 
erwise have very similar properties. 

We cannot exclude that other classes of AGN and 
QSOs emit bursts with a comparable rate since the derived 
3 a upper limits (see Table 6) are above the emission rate 
of the radio-quiet quasars (0.0146 y _1 per quasar). The 
lower emission rate of radio-loud quasars argues against a 
direct connection between radio jets and GRB emission. 

5.2. Comparison with previous results 

Two previous studies may appear in contradiction with 
our results. However, a thorough analysis shows that one 
of these is not comparable with our study, and that the 
other in fact does not contradict our results. 

5.2.1. The Kurt & Romberg (1995) analysis 

Kurt & Romberg (1995) compared the two-point angu- 
lar correlation function of 458 gamma-ray bursts with 
the two-point angular correlation function of a sample of 
413 unresolved extragalactic radio sources (QSO, BL Lac 
and radio galaxies). They found that the correlation func- 
tion of GRBs does not differ from a random distribution 
over all angular scales while the one for the radio sources 
shows significant autocorrelation on small angular scales 
at a level of > 3(7 . Kurt & Komberg conclude that GRBs 
do not belong to a normal population of galaxies with 
z < 1.5. This may not be correct in this generalization, 
particularly since these authors dealt with compact radio 
sources and did not consider radio-quiet objects. More- 
over, Drinkwater & Schmidt (1996) have recently pointed 
out the very different clustering properties of radio galax- 
ies and radio-loud quasars, both of which are well-mixed 
in Kurt & Komberg's sample. Because we found only co- 
incidences in excess of random expectation for radio-quiet 
quasars, our result is not in contradiction to that of Kurt 
& Komberg (1995). Furthermore, we found no hint for a 
correlation between radio-loud objects and GRBs which 
agrees with their conclusion. 

5.2.2. The Webber et al. (1995) correlation 

Webber et al. (1995) compared ~60 GRBs located to 
within very small error boxes of ~0.25 deg 2 or less with 
several catalogues of extragalactic objects. More than half 
of their sample of GRBs was taken from the first Interplan- 
etary Network (IPN) GRB Catalog (Atteia et al. 1987), 
seven GRBs from the COMPTEL catalogue of bursts 
(Hanlon et al. 1994) and 16 bursts from the third IPN 
(Hurley, private comm.). They find that any correlation 
between these catalogued extragalactic objects and the 
GRBs are consistent with chance expectation. In order 
to compare the bursts with the extragalactic objects they 
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count the number of the closest bursts falling in bins of 
width 0?5 in radial separation independent of the actual 
shape of the error box. They found two coincidences: QSO 
2219-420 (z = 1.3) in the error box of GRB 790419 and 
QSO 0116-288 (z = 0.798; which was discovered in a deep 
optical follow-up search of the GRB error box by Peder- 
sen et al. 1983) in the error box of GRB 781119 (obviously 
mis-printed as GRB 791119 in Webber et al. 1995). Web- 
ber et al. (1995) simulated 100 sets of quasar positions 
and obtained a background value of 3.7 ± 1.8. Motivated 
by the large increase of known quasars since the 1989 com- 
pilation of AGN and QSOs by Veron-Cetty & Veron as 
used by Webber et al. (1995), we reanalyzed Webber's 
sample of GRBs with the 1995 compilation of radio-quiet 
QSO. Because the error boxes of the seven bursts localized 
by COMPTEL have areas well within ~0.25 deg 2 , which 
is larger than most of the error boxes of the IPN GRB 
Catalog, we excluded the former here. For the 28 bursts 
from the first IPN GRB Catalog (Atteia et al. 1987), we 
compared the positions of radio-quiet QSO with the ex- 
act (elongated) shape of the error boxes and arrived at 
the same result as Webber et al. (1995): one quasar each 
is located in the error boxes of GRB 790419 and GRB 
781119. Then we determined the expected background 
value with 1000 simulations for the 28 burst positions. 
Accepting quasars within a bin width of 0?5 we derive the 
same result as reported by Webber et al. (1995), namely 
that a correlation between quasars and the GRBs is con- 
sistent with chance expectation. In a second approach, 
rather than using a standard radial distance of 0?5 for all 
bursts, we converted the area of a given error box to a cor- 
responding error circle for each burst. With this approach 
we expect 0.18 ± 0.4 random coincidences as compared to 
the above value of 3.7 ± 1.8. Thus, the probability to find 
two radio-quiet quasars in the error boxes of the GRBs 
is only 1.2%. This result corresponds to a probability of 
98.8% for a positive correlation, i.e. much higher than de- 
rived by Webber et al. (1995) and confirming our results. 

5.3. Host questions 

Our result of a correlation between GRBs and intrinsically 
bright, radio-quiet quasars does not necessarily imply that 
GRBs are produced in the nuclei of these quasars, but 
possibly in their host galaxies. Thus, there are two dif- 
ferent questions related to the hosts of GRB sources: Are 
host galaxies bright enough to be detectable in previous 
or future deep searches? What do we know about the host 
galaxies of radio-quiet quasars? 

5.3.1. Deep optical searches and the GRB hosts 

Earlier deep optical searches of small GRB error boxes 
(by imaging in bandpasses from U to K) did not find suit- 
able host galaxies (Schaefer 1992). This has become known 
as the "no-host" problem for cosmological GRB scenarios 



(Fenimore etal. 1993a). Recently, near-infrared imaging 
of six small GRB error boxes (a sample of GRBs different 
from that used by Schaefer (1992)) in the JHK' bands has 
shown that the "no-host" problem does not exist (Lar- 
son etal. 1996). On the contrary, even an overabundance 
of brighter (K < 15.5 mag) galaxies was found. Whether 
these might be hosts for radio-quiet quasars remains to be 
seen. 

As the detectability of quasars is concerned, the lack of 
excess quasars in GRB error boxes is neither contradicting 
our results nor surprising. In photometric studies, radio- 
quiet quasars are not separated from other types of AGN 
or radio-loud quasars, and thus the existence of any excess 
of radio-quiet quasars is easily missed. We should note 
here that a recent deep optical imaging (Vrba et al. 1995) 
has revealed marginal evidence for a quasar excess at the 
rate of one per high Galactic latitude GRB error box. 
Further studies able to distinguish between the different 
types of AGN are certainly needed to clarify this situation. 

In addition, there are two probable kinds of quasars 
which are difficult to detect in optical, photometric 
searches. First, there are only a few quasars showing 
Seyfert 2 type spectra. If the ratio of Syl/Sy2 type spec- 
trum quasars were equally large as the ratio of Syl/Sy2 
galaxies, then a large number of quasars will have been 
missed in optical studies. Ultra-luminous infrared galax- 
ies discovered by IRAS (Sanders et al. 1988) might be a 
second possibility to explain bursts lacking optical coun- 
terparts. These are suggested to be young quasars in a 
dust-enshrouded initial phase. 

5.3.2. The hosts of radio-quiet AGN 

Until quite recently radio-quiet quasars were thought to 
be associated mainly with spiral galaxies, and radio-loud 
ones with elliptical galaxies (see Bahcall et al. 1995 for an 
overview of references to previous, apparently contradict- 
ing results). However, very recent deep images taken with 
the Hubble Space Telescope suggest that quasars reside in 
a wide variety of environments, not confirming the above 
preconception. If radio-quiet quasars are indeed equally 
likely to be associated with both spirals and ellipticals, 
then their clustering properties are difficult to predict. To- 
gether with the reasonable assumption that GRBs are an 
extremely rare phenomenon in the life of a galaxy, the lack 
of autocorrelation among GRBs seems more natural. 

On the other hand, the production of GRBs must be 
related to some specific property of the quasar, otherwise 
one would expect also a correlation of GRBs with normal 
spirals, which is not observed. In the next section it is 
shown that the detected number of GRB and the total 
number of quasars is in agreement with the emission rate 
determined in section 4.2. Because the number of spirals 
exceeds the number of quasars by a factor of order 10 3 , 
our result argues strongly against a relation between GRB 
and normal spirals. 
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5. 5. Models for GRBs from quasars 
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In a Friedmann universe with zero cosmological constant 
the luminosity distance d; is defined as 

d L = {°o z + (<?o - 1) (1 + 2q Q zy - I 

whith Hq the Hubble constant, go the deceleration param- 
eter and c the velocity of light. The comoving volume V(z) 
can be approximated by: 

For go = 0.5 this correspond to: 



32?r 



c 

#0 



3 / [z + 1 - ' 




7.24 10 12 Mpc 6 



assuming Ho — 50 Mpc" 1 km s _1 . From Figure 2 in 
Hewctt ct al. (1993) the cumulative space density of 
quasars $(< Mb j; Mpc~ 3 ) for the optimized quasar sam- 
ple can be estimated. Assuming Mb,, ~ M v < —24.1 we 
get $ = 8 10~ 8 Mpc~ 3 for the redshift range 0.2 < z < 0.5 
and 2.2 10" 7 Mpc~ 3 for the redshift range 0.5 < z < 1.0, 
respectively. To derive the total number of quasars we 
further assumed that < &o.2<z<o.5 holds for redshifts below 
z = 0.2 and that 15% of the quasars are radio-loud. Under 
these assumptions we can estimate that in total 2.9 10 4 
quasars meet the criteria of the optimized quasar sample. 
The number of quasars studied, 967, therefore corresponds 
to 3.4% of the total sample. This number has to be com- 
pared with the percentage of GRBs cross-identified with a 
radio-quiet QSO. According to Table 8 we found counter- 
parts to 12.5 GRBs for a total sample of 134 bursts stud- 
ied, or a detection rate of 9.3%±2.8%. Thus, the expected 
detection rate is within two sigma of the observed one, a 
reasonable agreement given the small number statistics. 

One should also be aware that there are two biasing ef- 
fects which increase the rate of optical counterparts. The 
first is that more luminous quasars at a similar redshift 
show a higher probability to be detected in most kinds 
of surveys. Our finding that more luminous quasars have 
a higher probability of being a GRB counterpart suggest 
that the basic quasar sample used is already biased to- 
wards a higher rate of counterparts. A further bias towards 
a higher rate of counterparts is caused by the clustering 
of radio-quiet quasars (Shanks & Boyle 1994). Because 
the clustering occurs on small scales it does not affect our 
results given the large error boxes of the GRBs. On the 
other hand, unknown quasars show a higher probability 
to be located near to a known quasar in comparison to an 
arbitrary location on the sky. 



Merging neutron stars and black holes are one possibility 
to explain quasars as counterparts of GRBs and are dis- 
cussed in the first of the following subsections. There are 
not many publications explaining GRB with emission by 
the active nucleus itself. Because our result strongly sug- 
gests such a mechanism we show in a second subsection 
that the gamma-ray emission of radio-quiet quasars/ AGN 
can reach the energy range of GRBs. Finally, we discuss 
the model of Leiter (1980) which may explain the observed 
correlation. 



5.5.1. Merging neutron stars and black holes 

An interesting scenario to explain GRB from AGN and 
quasars was proposed by Epstein et al. (1993). Because 
compact-star binaries have a higher mass than normal 
stars they settle to the center of the nuclear stellar clus- 
ter (fraction of a parsec) before they coalesce. The binary 
merging will typically occur in a time span of a second 
and will accelerate debris to very high velocities reach- 
ing almost the speed of light. The quiescent emission of 
the quasar nucleus is then up-scattered to the range of 
gamma-rays by this accelerated hot material (debris), i.e. 
the spectrum is shifted by T 2 in photon energy and the 
flux is boosted by T 4 . The GRB spectra and their light 
curves are therefore explained by the relative orientation 
between the observer, the beam direction of the emission 
of the active nucleus, and the path of the accelerated ma- 
terial. This scenario predicts that the spectra and vari- 
ability of the GRBs would be correlated with that of the 
host AGN. 



5.5.2. Steady gamma-ray emission from quasars 

Already the HEAO-1 observations (Rothschild et al. 1983) 
revealed AGN to show emission up to the 100 keV energy 
range (mainly from Seyfert Is). Emission up to 100 MeV 
was detected from radio-loud quasars and BL Lac objects, 
for which synchrotron emission is the favoured mechanism. 
However, no emission above 500 keV was detected from 
radio-quiet AGN or QSOs with the three CGRO instru- 
ments BATSE, COMPTEL and OSSE (Schonfelder 1994). 
It appears that two different mechanisms are at work in 
the high-energy band depending on the radio-loudness of 
the objects. The existence of two mechanisms is supported 
by the well-established difference between radio-loud and 
radio-quiet QSOs seen in X-rays (cf e.g. Wilkes & Elvis 
1987; Schartel et al. 1996). Our result fits into this phe- 
nomenological picture, because we found GRB counter- 
parts only for radio-quiet quasars. On the other hand, the 
maximum energy of photons from GRBs reaches up to 
several GeV and thus is significantly higher than the ob- 
served steady 500 keV emission from radio-quiet AGN and 
QSOs. Thus, there is a clear difference in the gamma-ray 
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spectra of GRBs and the steady emission of radio-quiet 
quasars. 

Possibly, the MISO and Ariel 5 observations of NGC 
4151 (Pcrotti ct al. 1981), a radio-quiet Seyfert galaxy, 
may give some clues the understanding of our results. 
The flux above ~500keV from NGC 4151 was found to 
be variable by a factor of 3 to 10 on a time scale up to 
1 year as measured with Ariel 5 in December 1976, and 
MISO in May 1977 and September 1979. The observed 
7-ray fluxes are more than an order of magnitude higher 
than the limits reported from the OSSE/COMPTEL mon- 
itoring of NGC 4151 (Schonfeldcr 1994) or any previous 
measurement. Therefore, such a possible gamma-ray high 
state seems to be extremely rare. The variability by a fac- 
tor of 3 to 10 in the >500kcV range on a time scale of 
less than one year should be compared with a hard X- 
ray variability of only up to a factor of two (Perotti et al. 
1981). 

5.5.3. The model of Leiter (1980) 

One model to explain the variability found at the high- 
est measured energies was suggested by Leiter (1980). If 
a massive (M > 10 8 M Q ) Kerr black hole exists in the nu- 
cleus of NGC 4151, then the Penrose Compton-scattering 
mechanism can lead to the production of 7-rays with a 
cutoff at an energy of about 3 MeV. The burst is triggered 
by sporadic injection of hot plasma caused by turbulence 
in the inner region of the accretion disk. The resulting 
7-ray burst is emitted within an angle of 40° from the 
equatorial plane of the Kerr black hole. The duration of 
such an event is given by the size of the target region di- 
vided by the speed of light which is located between the 
event horizon and the ergosphere. In most unified models 
it is usually assumed that radio-loud quasars contain a 
rotating black hole matching the criteria of Leiter's sim- 
ulations. But for radio-quiet quasars black holes without 
rotation are generally favoured. However, the accretion 
of matter through an accretion disk permanently trans- 
fers angular momentum to the black hole. The increasing 
angular momentum causes an increase of the ergosphere. 
Finally, the ergosphere is larger than the event horizon 
creating a very small target region for Penrose Compton- 
scattering. Further investigations have strengthened the 
possibility to produce GRBs in hot tori around black holes 
(Jaroszynski 1996). 

If we assume that at least a part of the optical lu- 
minosity of the quasar is caused by turbulence in the 
accretion disk, then the positional coincidence we found 
between GRB and intrinsically bright quasars can be ex- 
plained with a higher injection rate of hot plasma and a 
higher transfer of angular momentum according to Leiter's 
model. Since particularly the X-ray emission of quasars is 
interpreted as a measure of the activity of the inner part of 
the accretion disk, a search for positional coincidences be- 
tween X-ray bright, radio-quiet quasars and GRBs would 



be most promising. It is interesting that the model pro- 
posed by Leiter already predicts a focusing orthogonal to 
the angular momentum of the black hole. This may ex- 
plain our negative result obtained for radio-loud quasars 
and BL Lac objects if we assume that the radio-jet is emit- 
ted along the rotation axis. 

6. Concluding remarks 

The aim of our study was to develop a method to handle 
the large error boxes of GRBs and to study the positional 
coincidences with QSOs and AGN. We found the sur- 
prising result that the nearest and intrinsically brightest 
radio-quiet quasars have a probability of >99.7% of being 
associated to the 134 best localized (and thus in general 
brightest and nearest) GRBs of the 3B catalogue. While 
previous claims of correlations have usually been found 
spurious in light of improved (and larger) data samples, it 
is not the fact of the mere correlation which is surprising, 
but rather the high probability we find. 

If our results are confirmed in future studies with im- 
proved statistics then bright quasars emit 7-ray bursts at 
an emission rate of 0.015±0.004 bursts per year and quasar 
for the BATSE sensitivity limit, based on a comparison of 
the 134 brightest BATSE bursts with 967 quasars of an 
"optimized quasar" sample. In order to determine the en- 
ergy emitted during a burst we consider only the 20 GRBs 
with a single quasar in their error box. For 16 of them the 
fluxes are given in the 3B catalog, thus yielding a mean 
total emission of 3.1 10 51 erg for the 20-50 keV band, of 
7.0 10 51 erg in the 50-100 keV band, and 3.2 10 52 erg in 
the 100-300 keV range, respectively. An isotropic emission 
over 4-7T sr was assumed and the following parameters were 
used: H o =50 km s _1 Mpc -1 , and q = 0.5 assuming power 
law spectra with a — 2.0. 

The deviation of the cumulative peak flux distribution 
of GRBs from the -3/2 slope should then be used with cau- 
tion to study the curvature of space-time, because quasars 
themselves show an intrinsic evolution. 

We neither wish to discuss the implications of our re- 
sults on models, nor to build a theoretical scenario of GRB 
emission from radio-quiet quasars. However, the brief re- 
view of the scenarios proposed by Epstein etal. (1993) 
and Leiter (1980) shows that there are already theoretical 
ideas matching our findings, and (in the case of Leiter's 
scenario) that phenomena related to neutron stars (or 
stellar-sized compact objects) are not the only possible ex- 
planation for GRBs. Also, the measurements of NGC 4151 
provide evidence for emission from radio-quiet AGN and 
QSO in the energy range required for GRBs. This implies 
that in the Epstein et al. (1993) scenario the boosting dur- 
ing the GRB is rather moderate. Thus, our finding of a 
positive correlation of GRBs with radio-quiet AGN com- 
bined with some previously suggested scenarios of cosmo- 
logical GRBs may provide a new avenue towards under- 
standing the GRB puzzle. 
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